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ABSTRACT 


This  report  presents  the  results  of  a  two-phase  research  program  consist¬ 
ing  of  (1)  in-flight  measurement  of  aerodynamic  pressures  and  structural 
loads  on  a  compound,  rigid-rctor  helicopter  and  (2)  correlation  of  these 
data  with  theoretical  results. 

Flight  test  data  obtained  in  Phase  I  and  recorded  on  an  oscillograph  were 
read  on  an  oscillograph  reading  machine  and  were  processed  in  an  automatic 
data  reduction  program.  This  data  processing  consisted  of  integration  of 
the  pressure  data  to  obtain  the  distribution  of  aerodynamic  lift  and 
pitching  moments  over  the  rotor  blade,  as  functions  of  azimuth  position. 
Airload  and  structural  load  data  were  harmonically  analyzed. 

Output  of  the  data  reduction  program  was  used  in  Phase  II  as  input  to  the 
correlation  program.  Hie  measured  airloads  were  used  to  compute  the  theo¬ 
retical  bending  and  torsion  responses  or  the  blade.  The  measured  torsion 
moments  were  used  in  the  theoretical  prediction  of  the  airloads.  Hie 
results  of  the  applied  theories  are  compared  with  the  flight  measurements. 
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FOREWORD 

This  report  describes  a  two-phase  research  program  consisting  of  (l)  flight 
test  measurements  of  helicopter  rotor  blarJe  structural  loads  and  aerody¬ 
namic  pressures  and  (2)  correlation  of  these  measurements  with  data  obtained 
from  current  theories.  This  research  program  was  conducted  by  the  Lockheed- 
California  Company  under  Contract  DA  Wt-i77-AMC-357(T)  to  the  U.S.  Army 
Aviation  Materiel  Laboratories  (USAAVLABS),  Fort  Eustis,  Virginia. 

The  research  program  was  performed  during  the  period  from  June  1966  to 
October  1967.  Technical  monitoring  of  the  project  for  USAAVLABS  was  by 
W.  E.  Hetties. 


The  report  covering  the  program  is  presented  in  three  volumes.  Volume  I 
is  entitled  "Measurement  and  Data  Reduction  c':  Airloads  and  Structural 
Loads".  It  contains  the  main  body  of  the  report  plus  Appendixes  I  through 
IV.  Volume  II  contains  Appendixes  V  through  IX,  with  all  flight  test  data 
in  tabular  form.  The  correlation  of  the  measured  airloads  and  structural 
loads  with  theoretical  data  is  covered  in  Volume  III,  "Theoretical 
Prediction  of  Airloads  and  Structural  Loads  and  Correlation  with  Flight 
Teat  Measurements". 

The  Lockheed  program  was  under  the  technical  direction  of  A.  W.  Turner  and 
W.  E.  Spreuer,  engineering  managers ,  and  J.  E.  Sweers,  project  leader.  The 
test  pilot  was  1!.  Goudey.  Additional  Lockheed  personnel  associated  with 
the  program  included  W.  H.  Foulke  and  R.  A.  Berry,  flight  test; 

C.  J.  Buzzetti,  E.  A.  Bartsch,  S.  H.  Lomax,  and  T.  H.  Oglesby,  structural 
flight  measurement ;  R.  H.  Cook  and  R.  G.  Murison,  instrumentation; 

R.  P.  Baker  and  V.  C.  Weddle,  data  processing;  R.  E.  Donham  and  D.  E.  Jan  da, 
rotary  wing  dynamics;  C.  H.  Ranschau,  programming;  and  R.  P.  Beal,  editor. 

Appreciation  is  due  USAAVLABS  for  their  help  in  providing  assistance  and 
advice  in  planning  and  executing  the  entire  research  program. 
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INTRODUCTION 


The  problem  of  analytically  predicting  loading  conditions  of  helicopter 
rotor  blades  in  high  speed  flight  consists  of  three  parts: 

•  To  find  the  conditions  for  trim,  i.e.,  vehicle  attitude  and  blade 
control  angles. 

•  To  find  the  distribution  of  airloads  over  the  blades  as  a  function 
of  azimuth  position. 

•  To  find  the  response  of  the  rotor  system  to  a  given  airload  dis¬ 
tribution. 

The  three  parts  of  the  analysis  are  interrelated,  since  the  trimmed  condi¬ 
tion  of  blade  control  angles  is  determined  by  the  airloads  developed,  and 
because  the  blade  hub  stiffness  is  normally  low  enough  to  result  in  blade 
motions  of  such  magnitude  as  to  affect  the  airloads  distribution  substan¬ 
tially.  Furthermore,  the  airloads  do  not  vary  linearly  with  the  control 
angles  or  the  blade  response.  Starting  with  a  given  weight  distribution 
and  flight  speed,  therefore,  requires  an  iteration  procedure  to  obtain  the 
structured  loads  on  the  blades. 

Since  in  each  of  the  steps  in  the  analysis  certain  approximations  and  simpli¬ 
fying  assumptions  must  be  made,  a  considerable  inaccuracy  in  the  final 
result  of  bending  and  torsion  moments  can  occur.  It  i .  therefore  of  great 
value  to  obtain  an  intermediate  result  for  comparison  such  as  measured  air¬ 
load  distributions. 

In  the  flight  tests  described  in  Volumes  I  and  II,  measurements  of  air  pres¬ 
sures  and  blade  stresses  are  used  to  determine  airload  distribution  and 
blade  response,  while  the  blade  control  angles  and  flight  attitude  are  also 
measured.  This  males  it  possible  to  separate  the  three  parts  of  the  analysis 
from  each  other  and  to  compare  the  results  of  each  separable  analysis  with 
test  data. 

In  the  theoretical  work,  extensive  use  was  made  of  the  Cornell  Aeronautical 
Laboratory's  computer  program  for  rotor  blade  loads  analysis;  in  addition, 
two  Lockheed  computer  programs  were  used. 
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ANALYTICAL  METHODS 


In  this  section  a  brief  description  is  given  of  the  theoretical  methods 
used  in  the  analysis  of  rotor  loads. 
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COMPUTER  PROGRAMS 


In  the  following  discussion  the  two  Lockheed  rotor  loads  analysis  programs 
used  in  the  theoretical  work  are  referred  to  as  Program  I  and  Program  II. 
Similarity  between  the  two  programs  makes  it  possible  to  give  a  parallel  de¬ 
scription  of  both. 

Program  I  performs  the  computation  of  rotor  blade  responses  from  a  given 
airload  distribution.  This  airload  distribution  consists  in  general  of  the 
harmonics  of  lift,  drag,  and  aerodynamic  pitching  moments ,  which  are  given 
as  lumped  loads  at  up  to  20  blade  stations.  Program  II  performs  similar 
coasputations  but  allows  the  use  of  a  larger  number  of  blade  stations.  In 
addition.  Program  II  includes  the  computation  of  aerodynamic  lift  and  drag 
at  these  blade  stations.  (A  flew  chart  of  Program  II  is  given  in  Figure  1.) 

In  the  airloads  computations,  the  following  assumptions  are  made: 

•  The  induced  velocity  distribution  can  be  obtained  from  a  descrip¬ 
tion  of  the  wake.  The  wake  is  approximated  by  a  number  of  Helmholtz 
ring  vortexes  representing  the  trailing  tip  vortexes  of  the  rotor 
blades.  The  vertical  displacements  and  the  diameters  of  the  ring 
vortexes  are  specified.  (The  vertical  velocity  of  propagation  of 

the  vortex  rings  and  the  contraction  depend  upon  wake-on-vake  effects. 
A  major  contribution  to  the  vertical  velocity  is  the  self-induced 
velocity  of  the  vortex  ring  which  is  a  function  of  vortex  strength 
and  core-radius  to  ring-radius  ratio  (Reference  2).  For  hovering 
conditions,  an  estimate  of  the  ring  spacing  may  be  obtained  from 
Reference  2.  At  low  forward  speed  these  estimates  are  assumed  still 
to  be  valid.  At  higher  forward  speeds  the  rotor  is  essentially 
unloaded  and  the  selection  of  ring  vortex  spacing  is  of  lesser 
importance.) 

•  The  strength  of  the  ring  vortexes  is  obtained  from  the  total  rotor 
thrust  and  is  taken  as  constant  around  the  azimuth.  The  shed  vortic- 
ity  is  neglected. 

•  The  induced  velocity  due  to  wake  vorticity  is  approximated  by  using  a 
number  of  simple  algebraic  expressions.  This  eliminates  the  very  time- 
consuming  evaluation  of  the  elliptic  integrals  in  the  expressions  for 
the  exact  velocities  resulting  from  the  Helmholtz  ring  vortex  model. 

•  Optionally,  the  induced  velocity  can  be  taken  as  uniform  over  the 
rotor  disc.  In  this  case,  it  is  obtained  from  momentum  theory. 


•  The  lift,  drag,  and  pitching  moment  coefficients  are  expressed  as 
functions  of  the  local  angle  of  attack  and  Mach  number,  as  veil  sis 
of  the  profile  thickness.  These  expressions  are  given  in  tabular 
form.  The  local  values  of  the  coefficients  are  found  by  interpola¬ 
tion. 

•  Initially,  the  rigid  blade  geometry  is  used.  The  effects  of  blade 
element  velocities  p*».r*0.1el  and  perpendicular  to  the  shaft,  due  to 
structural  deformations,  and  of  changes  in  local  coning  angle  and 
blade  twist  due  to  elasticity  are  introduced  in  the  next  cycle,  where, 
after  the  blade  responses  are  computed,  the  airloads  are  reevaluated. 
The  linear  aeroelastic  effects  sure,  however,  already  introduced  in 
closed  form  in  the  first  cycle  (i.e.,  without  iteration)  as  shown 
below.  (The  term  linear  is  used  here  to  indicate  that  the  aero- 
elastic  effects  acting  upon  each  harmonic  sure  produced  by  vibrations 
in  the  same  harmonic,  while  the  interaction  between  one  harmonic  and 
another  is  ignored  at  this  point. ) 

The  simplified  theory  for  induced  flow  due  to  ring  vortexes  is  obtained  as 
follows. 

The  potential  flow  induced  by  &  ring  vortex  of  infinitesimal  core  thickness 
can  be  solved  by  an  iterative  procedure.  The  three-dimensional  problem  is 
first  converted  to  a  two-dimensional  problem  by  I’iing  the  apparent  symmetry 
of  the  ring  vortex  model.  Then,  using  a  number  of  grid  points  in  the  plane  of 
symmetry ,  the  Laplace  equation,  expressing  zero  divergence,  is  solved  at  each 
of  these  grid  points  in  succession  (except  for  the  grid  points  in  the  plane 
of  the  ring  vortex,  where  the  potential  is  constant). 

The  flow  pattern  resulting  from  these  computations  can  be  visualized  sis  a 
system  oi  equipotential  lines  and  streanlines  which  intersect  each  other  at 
right  angles.  This  flow  pattern  is  unique,  i.e.,  it  is  independent  of  the 
vortex  strength  or  the  dimension  of  the  ring  vortex. 

In  order  to  arrive  at  a  set  of  simple  expressions  for  the  computation  of 
the  induced  flow  components  at  an  arbitrary  point  in  the  flow  field,  the 
flow  pattern  is  modified  by  approximating  the  equipotential  lines  as  well 
as  the  streamlines  by  circles,  still  intersecting  each  other  at  right  angles. 
The  equipotential  lines  are  centered  on  the  centerline  of  the  ring;  the 
streamlines  are  centered  in  the  plane  of  the  ring. 

From  the  geometry  of  this  simplified  flow  pattern,  the  induced  velocities  are 
easily  obtained  in  terms  of  an  arbitrary  constant.  This  arbitrary  constant 
is,  of  course,  a  linear  function  of  the  vortex  strength.  One  more  modifica¬ 
tion  is  needed.  It  appears  that  if  the  arbitrary  constant  is  adjusted  to 
give  the  correct  value  of  the  flow  velocity  in  the  vicinity  of  the  vortex 
filament,  the  flow  velocity  at  the  center  is  underestimated  oy  a  factor  n/2. 
The  constant  is  therefore  made  a  function  of  the  radial  distance  at  vhich  the 
streamline  passes  through  the  plane  of  the  vortex  ring.  This  function  is 
selected  in  such  a  manner  that  the  exact  flow  velocities,  both  at  the  center 
and  in  the  vicinity  of  the  filament,  are  obtained. 
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One  of  the  important  features  of  the  Lockheed  airloads  program  is  the 
possibility  of  trimming  the  rotor  to  a  given  set  of  integrated  airloads 
acting  on  the  shaft.  Up  to  three  forces  and  three  moments  can  be  selected  to 
which  the  rotor  can  be  trimmed.  These  are:  forward  force,  side  force,  rotor  j 

thrust,  roll  moment,  pitch  moment,  and  torque.  The  trimmed  conditions  are  ! 

obtained  by  adjusting  the  same  number  of  the  following  angles:  collective  1 

pitch,  cyclic  pitch  (two  components), shaft  angle,  flight  path,  and  yaw  j 

angle.  Theoretically  any  combination  of  loads  to  be  trimmed  to  and  angles  * 

to  be  computed  can  be  made,  with  the  obvious  restriction  that  their  respec¬ 
tive  numbers  must  be  the  same.  The  practical  limitations  of  the  method 
are  the  sane  as  those  of  the  physical  rotor  system,  e.g.,  in  hover,  changes 
in  the  shaft  angle  will  not  effect  any  changes  in  the  total  loads.  It 
should  also  be  noted  that  if  it  is  attempted  to  trim  the  re*or  to  a  load 
or  set  of  loads  beyond  its  physical  capability,  the  solution  of  the  computer 
program  will  not  converge,  and  the  computations  will  be  terminated. 

The  method  by  which  the  trimmed  conditions  are  obtained  is  briefly  described 
as  follows. 

In  addition  to  the  computation  of  the  airloads  on  the  blade  elements  using 
the  given  starting  values  of  the  angles ,  the  airloads  are  also  evaluated 
after  making  a  small  unit  change  in  each  of  these  angles  in  succession. 

The  total  rotor  loads  (3  moments  and  3  forces)  are  then  computed, and  those 
found  for  the  starting  values  are  subtracted  from  each  of  the  others.  This 
results  in  a  6  x  6  matrix  of  changes  in  total  rotor  loads  due  to  changes 
in  angles.  From  this  matrix  the  columns  pertaining  to  loads  to  which 
the  system  is  not  to  be  trimmed  are  removed.  Also  the  rows  pertain¬ 
ing  to  angles  which  are  not  used  for  trimming  (but  are  given  as  "fixed" 
quantities)  are  removed.  The  resulting  square  matrix  is  inverted  and 
transposed,  resulting  in  the  "trim  matrix". 

The  rotor  loads  computed  with  the  starting  values  of  the  angles  are  sub¬ 
tracted  from  the  rotor  loads,  specified  in  the  input,  to  which  the  rotor 
must  be  trimmed.  This  results  in  a  column  matrix  of  six  imbalances.  The 
imbalances  in  loads  to  which  no  trim  is  required  sure  removed. 

Finally,  the  column  of  imbalances  is  premultiplied  by  the  trim  matrix,  and 
the  resulting  column  is  scalar  multiplied  by  the  unit  of  change  of  the  angles, 
resulting  in  the  '■  inges  in  angles  required  for  trim.  These  changes  are 
based  on  linear  interpolation  or  extrapolation.  The  process  is  therefore 
repeated  until  convergence. 

After  convergence,  the  final  airloads  on  the  blade  are  computed .  These  do 
nc  .  yet  include  any  aeroelastic  effects.  The  linear  aeroelastic  effects 
a:  e  found  by  repeating  the  airloads  computations  after  making  successive 
unit  changes  in  vertical  (axial)  air  velocity,  tarjential  air  velocity,  j 

collective  pitch  angle,  collective  pitch  angular  velocity,  and  uniform  [ 

coning  angle.  The  difference  between  these  airloads  and  the  normal  air-  { 

loads  is  averaged  over  the  azimuth,  resulting  in  the  "hovering"  (linear)  j 

aeroelastic  effects  for  each  blade  station.  This  information  is  sub-  I 

sequently  used  in  the  response  calculations.  j 

| 
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From  this  point  on,  the  two  Lockheed  programs  are  similar,  and  the  follow¬ 
ing  description  applies  to  both.  However,  some  major  differences  exist 
and  will  be  indicated  where  they  occur. 

Before  entering  the  response  part  of  the  program,  certain  inertia  forces 
are  added  to  the  lumped  airloads  on  the  blade  elements.  These  are  the 
effects  of  vertical  acceleration,  pitch  and  roll  rate,  pitch  and  roll 
acceleration,  and  rate  of  change  of  rotor  rpm.  In  Program  II  the  linear 
accelerations  in  fore-aft  and  side  direction  are  also  included.  The  total 
lumped  loads  are  referred  to  as  external  loads  in  the  following  discussion. 

The  response  to  the  steady- state  part  of  the  external  loads  is  found  by 
an  iteration  process.  First,  the  bending  moments  about  two  axes  and  the 
torsion  moments  are  ccaputed,  using  rigid  blade  geometry.  For  the  angular 
position  of  the  structural  principal  axes  of  inertia  of  the  blade  section, 
the  combination  of  collective  pitch  angle  and  built-in  twist  is  used  (the 
cyclic  control  angle  is  ignored  at  this  point;  however,  the  effect  of  cyclic 
angle  is  linearized  and  used  later  as  it  affects  the  first  harmonic 
excitation).  The  deformed  shape  of  the  rotor  blade  is  then  computed  based 
on  the  computed  moments.  The  changes  in  deformation  result  in  changes  in 
centrifugal  forces  as  veil  as  in  a  new  geometry.  By  using  the  new  geometry 
and  centrifugal  •‘urces,  a  new  set  of  bending  and  torsion  moments  is  computed. 
This  process  is  repeated  until  convergence. 

Normally  the  above  process  would  not  converge  at  all,  but  rather  diverge 
very  rapidly.  Convergence  is  obtained  by  arbitrarily  reducing  the  amount 
of  change  of  the  moments  (in  particular  the  normal  bending  moment)  to  a 
fraction  of  the  computed  change  Since  at  convergence  the  changes  in 
moments  are  zero,  this  reduction  does  not  alter  the  final  results. 


The  response  to  oscillatory  airloads  (and  gyroscopic  forces  in  conditions 
with  pitch  and/or  roll  velocities)  is  obtained  from  the  Lagrange  equation: 


_d_/3T\  3T  3V  = 

dt  \3q /  3q  3q  * 


(1) 


where  T  is  the  kinetic  energy,  V  is  the  potential  energy  and  Q  is  the 
generalized  force  in  the  generalized  coordinate,  q.  Acceleration  forces 
are  obtained  from  the  first  term,  centrifugal  forces  are  obtained  from  the 
second  term,  and  Coriolis  forces  are  obtained  from  both  the  first  and  the 
second  term,  while  the  third  term  gives  rise  to  structural  forces. 


The  generalized  coordinates,  q,  used  in  the  analysis  are  primitive  bending 
and  torsion  modes  which  are  of  approximately  the  same  shape  as  the  natural 
modes  of  a  nonrotating  uniform  beam.  One  advantage  of  the  use  of  primitive 
mode  shapes  instead  of  natural  modes  is  that  the  natural  mode  shapes  are 
in  general  complex,  i.e.,  vertical  and  horizontal  displacements  are  not  in 
phase,  due  to  coupling  by  Coriolis  forces  between  these  displacements. 
Furthermore,  experience  shovs  that  differences  in  the  equilibrium  rxisition  of 
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the  blade  under  steady  load  can  result  in  significant  shifts  of  natural 
frequencies.  This  would  make  it  necessary  to  compute  new  natural  frequencies 
and  mode  shapes  for  each  condition,  which  can  be  avoided  by  working  with 
primitive  modes  only  (which  are  not  necessarily  orthogonal)  and  accounting 
for  all  the  coupling  terms. 


The  Lagrange  equations  are  expanded  into  matrix  form  by  first  writing  the 
displacements  in  terms  of  generalized  coordinates: 


3yi 


(a) 


then  by  writing  the  kinetic  and  potential  energy  in  the  rotor  system  in  terms 
of  displacements ,  velocities,  and  curvatures,  and  finally  by  taking  second 
derivatives  with  respect  to  q^  and 


as  I  ox-Lows: 


3  [~d  /3T  VI  _3_/ 3t\  3  /aV  \ 

!dt  W  /J"  as. \3%)  vv 


(3) 


to  find  the  (n,  m)  elements  of  the  matrices  A,  B,  C,  and  D  in  the  equation, 
(-u»2[AJ  +  iufl  [D]  -  n2  [B]  +  [C] )  (q)  =  {Q},  (fc) 


where  {Q}  is  the  column  of  generalized  forces. 


Feedback  of  rotor  response  into  aerodynamics  is  accomplished  as  follows. 

As  part  of  the  airloads  computations,  the  change  in  airload  due  to  a  change 
in  relative  velocities  and  angle  of  attack  is  included  at  each  radial  and 
azimuth  position.  These  changes  in  airload  are  harmonically  analyzed.  The 
".eroth  harmonic  part  is  used  to  define  the  linear  aerodynamic  damping  and 
stiffness  terms  in  the  above  matrix  equation,  and  these  terms  are  added  to 
matrices  D  and  C.  The  harmonics  of  the  airloads  due  to  blade  response 
cannot  be  used  in  a  closed- form  solution.  Therefore  the  complete  aerodynamic 
load  distribution  is  recomputed  using  the  blade  response  found  in  the  previous 
cycle. 

The  computation  of  the  total  bending  moments,  shear  loads,  and  torsion 
moments  at  each  azimuth  position  is  performed  as  follows: 

From  the  harmonics  of  displacements  in  the  primitive  modes,  the  harmonics 
of  displacements  in  x,  y,  and  z  directions  and  the  rotation  0  are  computed; 
for  example. 


where  i  is  the  primitive  mode  number  and  p  is  the  number  of  the  harmonic. 
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Hie  accelerations  x,  y ,  z,  and  9  are  found  from  x,  y ,  z,  and  0  by  multiplying 
by  -(p-1)2  n  2.  Note  that  p=l  indicates  the  steady  state,  p*2  indicates  the 
first  harmonic,  etc. 

Inertia  forces  are  computed  on  each  blade  element  from  the  accelerations. 
Similarly  the  velocities  x  and  y  and  the  resulting  Coriolis  forces  are 
computed.  The  inertia  and  Coriolis  forces  are  added  to  the  external  loads, 
resulting  in  total  forces  on  the  blade  elements. 

The  total  moments  are  computed  at  each  azimuth  position  from  the  total 
forces  on  the  blade  elements  and  the  instantaneous  blade  geometry.  Both 
the  forces  and  the  geometry  at  the  selected  azimuth  positions  are  obtained 
from  harmonic  synthesis. 

The  output  of  the  computer  program  consists  of  a  listing  of  the  two  bending 
moments,  torsion,  two  shears,  axial  load,  and  the  y,  z,  and  9  coordinates 
at  each  selected  station  and  azimuth  position.  In  addition  the  same  data 
are  given  in  harmonics  at  each  blade  station. 

Optionally  the  same  program  can  be  used  to  compute  natural  frequencies  and 
mode  shapes.  This  consists  of  the  solution  of  the  eigenvalue  problem: 

(-w2[A]  +  iuiQtD]  +  ICJ)  {q}  *  {0}  (6) 

If  the  linear  aeroelastic  effects  are  included  in  the  D  and  C  matrices, 
the  eigenvalues  u  consist  of  real  and  imaginary  parts.  The  modulus  of  u 
is  the  natural  frequency;  the  ratio  between  the  imaginary  part  and  the  real 
part  represents  the  damping  ratio. 

If  the  linear  aeroelastic  effects  are  not  included,  the  eigenvalues  are 
found  to  be  real. 

Hie  natural  mode  shapes  consist  of  x,  y,  and  z  displacements  and  rotation  9. 
If  the  blade  is  rotating*  the  y  and  z  displacements  are  coupled  through 
coriolis  forces  and  are  therefore,  in  general,  out  of  phase .  It  can  be 
shown  that  each  blade  element  describes  an  elliptical  path,  as  illustrated 
in  Figure  2. 

CORNELL  COMPUTER  PROGRAM 

This  program  i3  described  in  detail  in  Reference  3.  In  comparison  with  the 
Lockheed  method,  the  following  differences  are  of  interest: 

In  the  Cornell  program  the  geometric  blade  angles  are  treated  as  known 
quantities,  including  the  rotor  angle  of  attack.  It  is  therefore  not 
possible  to  compute  the  collective  and  cyclic  control  angles  and  the  rotor 
shaft  angle  required  for  a  given  trimmed  condition  in  a  single  computer 
ion.  (A  trimmed  condition  may  b*  defined  by  given  rotor  thrust,  shaft 
moments,  and  rotor  torque.) 


The  description  of  the  wake  of  the  rotor  blades  is  more  elaborate  in  the 
Cornell  program.  The  shed  vorticity  as  well  as  the  trailing  vorticity  is 
included  in  the  analysis.  At  short  distances  behind  each  blade,  the  trail¬ 
ing  vorticity  is  represented  by  a  number  of  distinct  (lumped)  trailing 
vortexes  rather  than  a  single  trailing  tip  vortex. 

The  wake  geometry  is  specified  by  input  values  of  the  wake  transport 
velocity,  which  can  be  selected  differently  at  each  station  and  azimuth 
position,  but  is  constant  with  time.  (The  azimuthal  differences  in 
wake  transport  velocities  are  restricted  by  the  i.etbod  in  which  this 
input  is  defined  to  five  harmonics  only.)  As  a  result,  the  wake  spirals 
in  this  program  do  not  contract  and  are  located  at  equal  distances  from 
each  other. 

Only  "flapping”  degrees  of  freedom  are  considered  in  the  Cornell  programs. 
Torsional  deformations  of  the  blade  are  specified  up  to  and  including  the 
fifth  harmonic  only. 

A  constant  lift  curve  slope  is  used,  which  is  cut  off  at  stall.  The  effect 
of  the  Kacb  number  is  not  included.  Drag  and  aerodynamic  pitching  moments 
are  not  computed . 

The  structural  responses  computed  in  the  Cornell  program  are  limited  to 
flapping  bending  moments.  These  are  found  as  the  summation  >f  the  bending 
moments  in  each  of  the  natural  modes. 
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APPLICATION  OF  THEORY 


COMPUTATION  OF  NATURAL  FREQUENCIES  AND  MOLE  SHAPES 

The  natural  frequencies  and  mode  shapes  of  the  instrumented  blade  are  mea¬ 
sured  by  a  shake  ter.t .  In  addition,  the  theoretical  mode  shapes  and  fre¬ 
quencies  are  computed,  using  Lockheed  Program  I,  for  the  cantilever- 
nonrotating  blade.  The  comparison  of  the  test  and  analysis  of  these 
frequencies  and  mode  shapes  is  mane  in  order  to  establish  the  validity 
of  the  mathematical  model  used  to  represent  the  rotor  system.  This  compari¬ 
son  is  shown  in  Figure  3. 

In  the  above  analysis  the  blade  is  subjected  to  Ig  gravity  forces  as  it  is 
in  the  shake  tests.  Ir»  order  to  find  the  frequencies  and  mode  shapes 
which  vill  be  used  in  the  Cornell  program,  the  computations  are  repeated 
without  the  lg  gravity  forces  and  with  the  normal  rotor  rpn.  The  coupled 
bending  mode  shapes  are  shown  in  Figure  2. 

AIRLOADS  COMP'TATICNS 


A  complete  list  of  the  flight  conditions  analyzed  is  contained  in  Table  I. 
Application  of  Lockheed  Computer  Program 

The  airloads  computations  performed  with  the  Lockheed  computer  program  are 
to  be  considered  separate  airloads  analyses,  i.e.,  the  airloads  are  com¬ 
puted  without  iteration  from  a  given  blade  geometry  including  the  struc¬ 
tural  responses  of  the  blades.  Sample  analyses  showed  that  torsional 
blade  responses  are  of  major  importance  in  the  determination  of  airloads 
on  the  blade  (Reference  4)  and  must  be  included  in  the  computations. 
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The  torsional  deformation  used  in  the  analyses  is  based  on  the  actual 
measurements  of  the  torsion  moment  on  the  blade  (at  a  radial  distance  from 
the  3haft  of  115  inches)  in  combination  with  the  effective  torsional  stiff¬ 
ness  of  the  blade  inboard  of  that  station.  The  resulting  twist  angles 
are  used  as  constants  over  the  blade .  For  the  purpose  of  the  separate 
airloads  analyses  the  computer  program  was  modified  in  order  to  accept 
the  input  of  the  torsion  moment  (at  station  115)  defined  at  36  azimuth 
positions.  Theoretically  this  provides  for  the  use  o.  all  the  harmonics 
of  this  torsion  moment  up  to  the  seventeenth.  However*,  it  must  be  realized 
that  the  distribution  of  the  torsion  moments  over  the  blade  span  in  the 
higher  harmonics  will  be  quite  different  from  that  in  the  lower  ones;  hence 
the  higher  harmonics  of  the  torsional  deformations  could  only  be  obtained 
correctly  if  different  effective  stiffnesses  were  used  in  connection  with 
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VERTICAL  DEFLECTION  -  INCHES 


TABLE  I.  FLIGHT  TEST  CONDITIONS 
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the  higher  harmonic  torsion  moment? .  Furthermore,  in  the  higher  harmonics 
the  angular  displacements  could  be  expected  to  be  different  at  the  dif¬ 
ferent  stations. 

The  harmonics  of  the  torsion  moments  at  station  115  whicn  are  used  in  the 
separated  airloads  analysis  are  shown  in  Volume  II.  Ccmnarison  of 
the  harmonic  components  at  station  115  and  station  135  confirms  the 
above  observations.  The  absence  of  any  clear  phase  relationship  between 
the  torsion  moments  at  these  stations  indicates  that  more  than  one  mode 
participates  in  the  torsion  moment  response.  These  uncertainties  with 
regard  to  the  finer  details  of  the  distribution  of  the  torsional  deforma¬ 
tions  must  be  borne  in  mind  in  the  review  of  the  comparisons  of  the  results 
of  the  separated  airloads  analyses  with  the  measured  data. 

No  provisions  were  made  to  introduce  a  given  normal  or  in-plane  bending 
response  into  the  separated  airloads  analysis.  That  the  effect  of  in-plane 
bending  on  the  airloads  can  safely  be  neglected  is  obvious.  The  effect 
of  the  normal  bending  response  is  evaluated  on  a  sample  basis  by  two 
different  methods: 

•  In  addition  to  the  predicted  airloads,  the  Lockheed  computer 
program  provides  the  change  in  lift,  drag,  and  pitching  moment 
due  to  a  unit  incremental  relative  vertical  velocity,  ^n  estimate 
of  the  incremental  relative  vertical  velocity  due  to  blade  bend¬ 
ing  response  can  be  obtained  from  the  response  analyses . 

•  The  Cornell  computer  program  does  include  the  effect  of  flapping 
bending  on  the  airloads.  Comparison  of  airloads  computed 

with  this  program  using  three  flapping  modes  and  a  modification 
where  the  displacements  in  the  natural  modes  are  zeroed  /it  shows 
that  the  effect  of  blade  bending  response  upon  the  airloads  is 
relatively  small. 

I**  the  Lockheed  computer  program  the  effect  of  not  including  the  blade  bend¬ 
ing  in  the  separated  airloads  analysis  is  largely  compensated  for,  as  far  as 
the  first  harmonic  is  concerned,  by  making  slight  changes  in  the  cyclic 
control  angles,  as  will  be  discussed  presently. 

Since  the  Lockheed  programs  prov':  de  an  option  to  trim  the  rotor  to  various 
selected  forces  and  moments,  it  is  interesting  to  look  at  the  advantages  of 
using  this  option. 

Of  the  angles  measured  in  the  flight  test,  the  collective  and  cyclic  control 
angles  are  known  with  sufficient  accuracy;  the  rotor  .^gle  of  attack, 
however,  is  not,  for  a  number  of  reasons.  Among  these  are  the  possible  effect 
of  the  wing  and  body  on  the  free-stream  velocity  at  the  rotor  disc  and  the 
fact  that  bending  displacements  of  the  rotor  blade  elements  are  neglected 
(the  first  harmonic  of  these  changes  the  tip  path  plane). 
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Small  variations  in  the  cyclic  control  angle  can  compensate  for  neglecting 
the  first  harmonic  of  the  vertical  velocities  of  tne  blade  elements  due 
to  blade  bending. 

The  total  rotor  thrust  is  found  from  integration  of  the  airloads  on  the 
rotor.  The  pitch  and  roll  moments  on  the  rotrr  can  also  be  obtained  in 
this  way.  It  has  been  found,  however,  that  ■‘..cse  moments  do  not  agree 
very  well  with  the  moments  derived  from  the  first  harmonic  of  normal  bend¬ 
ing  at  station  6  in  combination  with  gyroscopic  and  inertia  moments  based 
on  measured  roll  and  pitch  rates  and  accelerations. 

It  may  be  assumed  that  the  measurements  of  the  normal  bending  moments  at 
station  6  are  highly  reliable.  Furthermore,  substantial  errors  in  the 
measurements  of  pitch  and  roll  rates  and  pitch  and  roll  accelerations  are 
not  to  be  expected.  Therefore,  in  view  of  the  fact  that  a  very  small  error 
in  the  measured  and  integrated  airload  distribution  can  result  in  a  sub¬ 
stantial  error  in  computed  shaft  moments,  the  trimmed  solutions  for  the 
airloads  computations  are  based  upon  the  measured  first  harmonic  of  normal 
bending  at  station  6,  extrapolated  to  the  shaft,  multiplied  by  2  to  account 
for  four  blades,  and  combined  with  the  inertia  and  gyroscopic  moments  on  the 
rotor  due  to  the  measured  pitch  and  roll  rates  and  accelerations . 

The  situation  resulting  from  this  method  of  analysis  is  that,  as  far  as 
the  steady-state  and  first  harmonic  airloads  are  concerned,  the  analytical 
results  are  believed  to  be  quite  accurate.  The  comparison  of  test 
and  theory  must  therefore  be  considered  as  a  check  on  the  performance  of 
the  airloads  measurements. 

In  addition  to  the  trimmed  analyses  described  above,  analyses  were  made  in 
which  the  rotor  was  not  trimmed.  In  the  untrimmed  rotor  analyses,  the 
measured  control  angles  and  anule  of  attack  were  used,  Mso,  the  effect 
of  the  blade  torsion  was  left  out. 

The  Lockheed  program  provides  for  input  of  selected  vertical  displacements 
and  contraction  of  the  vortex  rings .  In  the  untrimmed  analyses  an  assumed 
relationship  between  downward  displacement  and  wake  contraction  is  used, 
based  on  experimental  work  (Reference  5).  The  downward  displacement  is 
found  from  the  uniform  induced  velocity  based  on  momentum  theory.  This  then 
determines  the  contraction,  i.e.,  the  wake  radius.  To  the  downward  displace¬ 
ment  is  then  added  the  vertical  component  of  the  free-stream  velocity 
based  on  the  measured  angle  of  attack. 

In  the  trimmed  analyses  the  vertical  displacement  Z  (positive  down)  and  the 
vertex  radii  £  were  selected  as  follows: 

Conditions  Z  =  5,  12,  20,  28,  36,  UU 

1,  !*,  19  £  =  205,  200,  195,  190,  185,  180 

Conditions  Z  -  8,  25,  1*0 

5,  8,  11  £  =  205,  200,  195 
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Condition  l6 


Z  =  -10,  -20,  -30 

^  =  205,  200,  195 

Conditions 

21,  23,  25,  26,  27,  31,  33,  Z  =  U,  10,  20 

36,  37,  39,  1*0,  U6,  50  £  =  205,  200,  195 

Application  of  Cornell  Program 

In  the  analyses  made  using  the  Cornell  A.L.  computer  program  for  rotor  loads, 
the  following  selections  were  made  regarding  the  wake  description,  the 
torsional  response,  and  the  blade  dynamic  properties  and  geometry: 

A  very  important  element  in  the  computation  of  rotor  loads  using  wake  vorticity 
is  the  placement  of  the  wake  elements  with  respect  to  the  rotor.  In  the 
Cornell  program,  the  location  of  the  wake  elements  is  obtained  from  selected 
inputs  describing  the  downward  velocity  of  the  wake  elements.  In  all  com¬ 
putations  made  with  the  Cornell  program,  this  velocity  was  taken  as  the  uni¬ 
form  inflow  velocity  obtained  from  momentum  theory  combined  with  the  component 
of  the  free-stream  velocity  perpendicular  to  the  plane  of  the  rotor.  The 
latter  is  found  using  the  measured  angle  of  attack.  The  wake  is  further 


described  as  follows: 

•  Number  of  azimuth  segments  of  the  wake  mesh 

behind  each  blade  (i.e.,  number  of  shed  vortices)  =  2 

•  Number  of  trailing  vortexes  in  the  rolled-up 

wake  (tip  vortex  only)  =  1 

•  Number  of  wake  revolutions  =  3 

•  Wake  advance  (i.e.,  fraction  of  A?  =  2*/no.  of 

blades,  by  which  the  wake  is  advanced  azimuthally)  =0.7 

•  Distance  rolled-up  tip  vortex  is  moved 

inboard  from  the  blade  tip  =  30  inches 


(Note:  The  radius  of  the  tip  vortex  is  then  180  inches.) 

The  torsional  response  is  specified  in  the  input.  The  Cornell  program  per¬ 
mits  a  separate  (different)  input  for  the  harmonics  of  the  geometric  blade 
angle  at  each  blade  station.  The  responses  used  in  the  Cornell  program  were 
obtained  from  Lockheed  Program  II  and  were  computed  from  the  measured  airloads. 
Since  provisions  are  made  for  the  input  of  the  first  five  harmonics  only,  the 
higher  harmonics  of  the  torsional  response  are  not  included. 

The  cone  angle  of  the  XH-51A  rotor  blade  is  not  constant  but  varies  between 
3.2  degrees  on  the  inboard  blade  and  2.2  degrees  on  the  outboard  blade.  In 
the  Cornell  program  a  single  cone  angle  must  be  selected.  This  cone  angle 
was  taken  as  the  average,  2.7  degrees. 


16 


Since  the  rotor  angle  of  attack  is  not  found  in  the  test  data  with  sufficient 
reliability,  and  since  this  angle  cannot  be  found  from  an  internal  trim  pro¬ 
cedure,  the  trimmed  angle  of  attack  computed  with  Lockheed  Program  II  is  used. 

The  natural  frequencies  and  mode  shapes  computed  with  Lockheed  Program  1  are 
used. 

RESPONSE  ANALYSES 


Three  different  methods  are  used  to  compute  the  blade  responses  to  the 
given  measured  airloads.  It  may  be  noted  that  since  the  total  airloads 
are  known,  the  damping  forces  and  aeroelastic  effects  are  included  in  the 
excitation.  Therefore  the  only  damping  used  in  the  computations  is 
structural  damping,  which  is  estimated  at  .02  [complete  damping,  i.e.,  the 
stiffness  matrix  C  is  multiplied  by  (i  =  .02i)j. 

The  response  analyses  are  again  considered  as  "separated  analyses"  since 
no  feedback  to  airloads  of  the  response  need  be  considered. 

The  computations  are  performed  with  Lockheed  Program  I  ar  described  in  the 
previous  chapter.  In  these  computations  7  harmonics  are  used  and  18  azi¬ 
muth  positions  (20*  increments  in  azimuth).  Thirteen  primitive  modes  are 
used,  consisting  of  five  normal  bending,  five  in-plane  bending,  and  three 
torsion  modes. 

The  same  computations  are  performed  with  a  simplified  program.  The  difference 
between  t'.is  new  method  and  the  method  above  is  that  the  responses  in  the 
modes  are  computed  for  all  harmonics  including  the  steady  state,  thus  ignoring 
the  nonlinearities  in  the  steady-state  response.  The  normal  bending  moments 
are  now  found  from  a  summation  of  bending  moments  in  each  of  the  modes.  In 
these  computations  the  same  13  primitive  mode  shapes  are  used  as  above. 

Also  as  above,  the  equilibrium  position  of  the  blade  is  taken  as  the  stead** 
state  deformed  shape  under  the  steady-state  part  of  the  load. 

Since  it  is  not  practical  to  modify  the  Cornell  program  in  such  a  way  that  a 
given  airload  distribution  can  be  given  as  an  input, the  Cornell  method  for 
the  computation  of  normal  bending  moments  is  simulated  by  a  variation  of  the 
modified  program.  This  consists  of  using  only  five  primitive  modes  in  normal 
bending  and  of  using  the  undeflected  blade  shape  as  equilibrium  position. 

The  measured  distributions  of  differential  pressures  do  not  provide  any  in¬ 
formation  with  regard  to  the  chord  forces.  The  integrated  pressures  (over 
the  chord)  give  the  normal  force  and  aerodynamic  moments  only. 

In  view  of  the  large  chordwise  bending  moments  measured,  in  particular  at 
the  blade  root,  it  is  apparent  that  the  aerodynamic  chord  forces  cannot  be 
neglected.  Furthermore,  coupling  between  in-plane  and  normal  bending  results 
in  an  effect  of  the  chordwise  airloads  on  the  normal  bending  as  well  as  on  the 
in-plane  bending. 
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For  the  purpose  of  the  present  investigation,  a  modification  was  made  to 
Lockheed  Program  I  by  which  the  chordvise  airloads  are  derived  from  the 
normal  forces  as  follows: 

•  The  normal  force  coefficient  is  obtained  from  the  measured  normal 
force  and  the  computed  relative  wind  at  each  blade  section. 

•  The  chord  force  coefficient  is  found  as  a  function  of  normal  force 
coefficient  and  Mach  number.  This  relationship  is  shown  in  Figure  h. 

•  The  chordwise  airload  on  each  blade  element  is  then  computed  from 
the  relative  wind  velocity  and  the  chord  force  coefficient. 

The  lumped  normal  forces  and  chord  forces  at  each  blade  station  and 
azimuth  position  are  resolved  in  the  directions  parallel  and  perpendicular 
to  the  shaft  axis.  The  angle  of  rotation  used  here  is  the  instantaneous 
angle  found  from  the  collective  and  cyclic  control  angles 

8  =  0  +0,  cos  'f  +  0,  sin  ? 

o  1c  Is 


HARMONICS  OF  ROTOR  LIFT 


In  the  computations  of  the  responses  to  measured  airloads,  using  Lockheed 
Program  I,  the  harmonics  of  normal  shear  force  at  each  blade  station  are 
computed.  It  is  assured  here  that  the  responses  of  the  noninstrumentod  blades 
are  the  same  as  those  of  the  instrumented  blade.  From  this  it  follows  that 
the  total  •vertical  shear  transmitted  to  the  shaft  contains  only  the  fourth, 
eighth,  twelfth,  etc.,  harmonics,  since  these  are  the  only  harmonics  to  which 
the  blades  respond  in  phase  with  each  other. 

Of  the  sever  harmonics  used  in  the  response  analysis,  only  the  fourth  is 
therefore  considered.  In  Table  II  the  cosine  and  sine  components  and  the 
amplitude  of  the  fourth  harmonic  of  rotor  lift  are  given.  These  forces  con¬ 
sist  of  the  summation  of  the  fourth  harmonics  of  airloads  and  inertia  forces 
on  the  blade  element"  •.  multiplied  by  four  to  account  for  the  four  blades . 

The  cosine  component  is  the  instantaneous  force  acting  at  the  time  when  the 
blades  are  ato*,  90*,  180*  and  270*  azimuth  positions;  the  sine  component  is 
the  instantaneous  force  acting  at  the  time  when  the  blades  are  at  22.5*,  112.5* 
202.5*  and  292.5*  azimuth  positions. 
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CHORD  FORCE  COEFFICIENT 


Figure  Chord  Force  Coefficient  Versus  Normal  Force  Coefficient, 
0012  Profile 
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TABLE  II.  FOURTH  HARMONIC  COMPONENTS  OF  ROTOR  LIFT 


INCLUDING  ROTOR  INERTIA  FORCES 


Cosine 

Component 

Sine 

Component 

Ampl itude 

Condition 

Cosine 

Component 

Sine 

Component 

Amplitude 

26 

-120 

-16 

120 

-6U 

-128 

11*1* 

27 

-332 

-536 

632 

64 

1* 

61* 

28 

-  60 

-  28 

68 

0 

6o 

6o 

29 

-168 

-200 

260 

-32 

-108 

112 

30 

-316 

-368 

488 

-  8 

-  96 

96 

31 

-400 

-740 

840 

8 

-152 

152 

224 

-  52 

228 

33 

-24 

-  56 

64 

-108 

-321* 

3U1* 

31* 

-56 

-104 

116 

61* 

-261* 

272 

35 

-200 

-400 

448 

-12 

-288 

288 

36 

-112 

-128 

168 

28 

-296 

300 

37 

-  88 

-124 

152 

1*8 

-272 

276 

38 

-  84 

-120 

148 

140 

-2l*0 

260 

39 

-  56 

-  8 

60 

172 

-292 

3l*0 

1*0 

-  12 

-  96 

96 

-12 

76 

76 

41 

-  84 

-  84 

120 

-80 

-101* 

132 

42 

-100 

-104 

144 

-76 

-  68 

101* 

43 

-  44 

-120 

128 

-12 

-  81* 

81* 

44 

-172 

-  76 

188 

96 

-  1*0 

10l* 

45 

-  48 

-  96 

108 

-12 

-  6o 

60 

46 

-  36 

-  64 

76 

1*12 

-52l* 

668 

47 

-124 

-  32 

128 

-32 

12 

32 

48 

-  48 

-104 

112 

-60 

-  96 

116 

49 

-  68 

-  56 

88 

-61* 

-  32 

72 

50 

-  76 

-  24 

80 

l 
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COMPARISON  OF  THEORETICAL  RESULTS  WITH  TEST  DATA 

Comparison  of  measured  and  computed  responses  (Figures  5  through  102)  shows 
that  on  the  outboard  part  of  the  Made  the  harmonics  of  the  bending  moments 
can  be  obtained  with  reasonable  accuracy.  Notable  differences  between  test 
and  theory  exist,  however,  in  the  following  areas: 

•  At  the  inboard  stations  the  normal  bending  moments  computed  with 
Lockheed  Program  I  show  in  many  conditions  a  substantial  difference 
in  first  harmonic  content.  This  difference  may  be  attributed  to 
inaccuracy  of  the  airloads  measurements. 

•  some  conditions  the  steady-state  part  of  the  normal  bending 
moments  shows  a  considerable  difference  between  test  and  theory. 

•  In  several  conditions  the  agreement  between  measured  and  computed 
in-plane  bending  moments  is  very  poor.  This  indicates  th»t  the 
relationship  between  normal  force  and  chord  force  coefficients, 
used  to  compute  the  drag  loads  (Figure  k)  t  may  be  deficient  in 
certain  areas.  This  is  in  particular  the  case  for  the  larger 
angles  of  attack;  as  shown  in  Figure  6,  condition  U,  which  is  a 
hovering  pullup,  the  drag  loads  appear  to  be  overestimated. 

•  The  agreement  between  measured  ar.d  computed  torsion  moments  is 
rather  poor,  particularly  at  the  lower  flight  speeds.  The  computed 
torsion  moments  show  a  large  fourth  harmonic  response.  This  is  to 

be  expected,  based  on  measured  as  well  as  computed  torsion  frequencies. 
It  appears  that  the  mechanism  by  which  the  torsion  moments  are 
produced  is  not  fully  known  at  this  time.  Improvements  in  the 
analytical  methods  should  include  the  degrees  of  freedom  of  the 
control  system  and  a  more  refined  definition  of  the  blade  elastic 
axis  and  local  eg. 

•  Comparison  of  results  from  Lockheed  Program  I  and  the  modified 
Program  I  shows  the  improvement  which  is  obtained  by  computing  the 
bending  moments  from  the  summation  of  airloads  and  inertia  forces, 
rather  than  from  the  summation  of  modal  responses. 

Comparison  of  measured  and  computed  airloads  shows  that  in  most  cases  the  basic 
shape  of  the  airload  distribution  can  be  predicted  (Figures  103  through  176). 

For  the  20  conditions  for  which  a  full  data  analysis  is  available,  the 
trimmed  rotor  loads  are  computed  using  Lockheed  Program  II.  The  shaft 
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moments  in  these  conditions  are  obtained  from  extrapolation  of  the  first 
harmonic  of  normal  bending  at  station  6,  and  include  the  gyroscopic  and 
'■nertia  moments  on  the  rotor  computed  from  the  measured  pitch  and  roll  rates 
and  accelerations.  The  rotor  lift  to  which  the  rotor  is  trimmed  is  taken  as 
the  integrated  airload.  The  angles  used  in  the  trimming  procedure  are  the 
cyclic  control  angles  and  the  rotor  angle  of  attack,  except  in  conditions  1, 

4,  and  19,  which  are  essentially  hovering  conditions. 

The  rotor  angle  of  attack  found  in  the  trimming  procedure  with  Lockheed 
Program  II  is  used  in  the  same  20  conditions  in  the  Cornell  program  (Figures 
123  through  142).  The  other  analyses  with  the  Cornell  program  are  performed 
using  the  angle  o*“  attack  as  measured  in  the  flight  tests  (Figures  143  thro.ig'-. 
171). 

In  addition,  the  conditions  for  which  a  full  analysis  is  made  are  analyzed 
untrimmed,  using  Lockheed  computer  Program  II  (Figures  105  through  122). 

Comparison  of  the  trimmed  and  untrimmed  conditions  shows  the  usefulness  of  the 
trim  procedure.  This  also  indicates  that  the  measured  angle  of  attack  cannot 
be  used  for  the  computation  of  the  rotor  airloads.  In  some  conditions  this  is 
quite  obvious.  In  condition  l6  for  instance  (Figure  108),  which  is  an  auto¬ 
rotation,  the  flight  path  angle  is  definitely  nonzero  and  should  be  included 
in  the  analysis.  The  trimmed  solution  for  this  condition  resulted  in  a  rotor 
angle  of  attack  of  14.5®  relative  to  the  free-stream  velocity.  It  is  also 
interesting  to  note  that  the  computed  rotor  torque  in  the  triraned  condition 
turned  out  tc.  be  a  small  negative  value,  as  is  to  be  expected  in  autorotation. 

In  condition  19  (transition  to  forward  flight),  an  unsucessful  attempt  was 
made  to  obtain  a  converging  solution  by  trimming  to  the  measured  rotor  thrust 
using  a  variation  of  the  rotor  angle  of  attack  at  a  forward  speed  of  20  knots. 
The  exact  condition  is  not  known,  mainly  because  the  forward  speed  was  not 
specified.  A  converging  solution  was  obtained  by  using  the  collective  control 
angle  for  trim.  The  rotor  angle  of  attack  was  estimated  as  -20®,  and  the 
forward  speed  was  estimated  as  40  knots.  The  airloads  are  shown  in  Figure  109 
The  trimed  collective  control  angle  was  found  to  be  11.13°  as  compared  to  the 
measured  angle  of  11.09°. 

A  good  agreement  was  obtained  using  Lockheed  Program  II  for  the  hovering 
condition,  as  shown  in  Figure  103.  Figure  104  shows  the  effect  of  introducing 
a  small  forward  velocity  upon  the  airloads.  This  forward  velocity  was  taken 
as  10  knets,  which  appears  to  be  too  high  since  it  results  in  computed  loads 
in  the  forward  quadrants  at  the  outboard  stations  which  are  too  high.  However 
at  the  in-board  stations  a  reasonable  agreement  with  the  test  was  obtained. 

Encouraging  results  were  obtained  with  the  Lockheed  Program  II  in  trimmed 
analyses  for  the  lower  velocities  (conditions  5,  3;  see  Figures  105  and  106). 
The  ract  that  these  analyses  showed  better  agreement  with  the  tests  than 


those  with  tr.e  Corner program  is  attributed  to  an  improvement  In  the 
description  of  the  wake  geometry  (see  Application  of  Theory,  Airloads 
Computations ,  for  tne  values  of  vertical  displacements  and  radii  of  the  ring 
vortexes ) . 

All  flight  conditions  of  Table  I  were  analyzed  with  the  Cornell  A.L.  program 
of  Reference  1. 

Comparison  of  the  tests  and  Cornell  theory  shows  that  at  least  at  the  outboard 
blade  stations  a  good  agreement  is  obtained  in  the  shape  of  the  airload  curves 
at  moderately  high  forward  speeds  (see  Figures  131,  137,  138,  and  lUl). 

However,  at  the  low  speeds  this  agreement  is  very  poor,  as  showr  in  Figures  125 
and  127.  It  is  likely  that  the  computed  deformation  of  the  wai.e  (as  determined 
by  the  vertical  wake  transport  velocity)  should  not  be  based  on  the  uniform 
inflow  velocity  computed  from  momentum  theory  at  the  lower  speeds.  Particu¬ 
larly  in  hovering  conditions  (conditions  1  through  U),  the  computed  steady-state 
part  of  the  loads  at  the  inboard  stations  is  not  in  agreement  with  the  test 
measurements . 

In  the  course  of  the  analytical  work  it  became  increasingly  difficult  to  obtain 
convergence  of  the  Cornell  program  at  the  higher  forward  speeds.  Ho  con¬ 
vergence  was  obtained  in  conditions  27,  29.  30,  and  31. 

Finally,  for  conditions  5,  8,  25,  31,  and  37,  the  airloads  were  computed  in 
trimmed  flight,  without  an  external  input  of  torsional  responses.  In  these 
analyses  the  responses  were  computed  from  the  computed  airloads  and  ted  back 
into  the  airloads  computations.  The  results  are  shown  in  Figures  172  through 
176. 


CONCLUSIONS 


1.  In  most  of  the  conditions  analyzed,  the  computed  bending  moments  were 
found  to  be  in  good  agreement  with  the  meajursd  moments.  The  agreement 
between  computed  and  measured  torsion  moments,  however,  was  poor,  indicating 
that  improvements  are  required  in  the  mathematical  model  as  far  as  uhe 
torsional  properties  are  concerned.  (These  may  include  the  introduction  of 
the  control  system  degrees  of  freedom  in  the  model . ) 

2.  Comparison  of  measurements  with  airloads  obtained  with  the  Cornell  program 
shows  good  agreement  of  the  variations  of  the  airloads  over  the  azimuth. 

'-3.  At  the  lower  forward  speeds  the  rotor  wake  can  be  represented  by  a  set  of 
ring  vortexes. 

4.  \n  urgent  need  exists  for  a  method  for  predicting  wake  deformation  and  for 
programming  this  method. 

•5.  Torsional  responses  are  important  in  the  computation  of  airloads. 

J&.  A  rotor  trim  procedure  included  in  a  rotor  loads  program  has  been  shown 
to  be  useful  in  the  prediction  of  rotor  loads. 

7 •  Based  on  the  results  obtained  for  chordwise  bending  moments,  the 
relationship  between  normal  force  coefficient  and  chord  force  coefficient  used 
in  Lockheed  Program  I  should  be  reexamined,  particularly  in  high  angle-of- 
attack  areas.  Such  investigation  may  result  in  improved  estimates  of  drag 
coefficients. 
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IMAL  BENDING 


NORMAL  BENDING  MOMENT  -  IN. -LB  x  1 


Figure  9.  Response  to  Measured  Airloads  (13  Modes), 
Condition  11  (Left  Turn  At  8U  Knots  TAS) 
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Figure  11.  Response  to  Measured  Airloads  (13  Modes), 
Condition  16  (Autorotation) 
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Figure  18.  Response  to  Measured  Airloads  (13  Modes) 
Condition  31  (Level  Flight  at  232  Knots  TA3) 
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Figure  22.  Response  to  Measured  Airloads  (13  Modes ), 
Condition  39  (Pullup  at  206  Knots  TAS) 
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Figure  2k .  Response  to  Measured  Airloads  (13  Modes ) » 
Condition  1*3  (Pullup  at  163  Knots  TAS) 
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Figure  33.  Response  to  Measured  Airloads  (13  Modes), 
Condition  9  (Left  Turn  at  Jl  Knots  TAS) 


Figure  34.  Response  to  Measured  Airloads  (13  Modes), 
Condition  10  (Right  Turn  at  38  Knots  TAS) 
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NJRMAL  BENDING  MOMENT  -  IN. -LB  X  1<T  NORMAL  BENDING  MOMENT  -  IN. -LB  x  1 


Figure  35.  Response  to  Measured  Airloads  (13  Modes), 
Condition  13  (Collective  Pullup  at  84.5  Knots  TAS) 
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figure  36.  Response  to  Measured  Airloads  (13  Modes), 
Condition  l4  (Collective  Pullup  at  86  Knots  TAS) 
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Figure  37*  Response  to  Measured  Airloads  (13  Modes), 
Condition  15  (Collective  Pullup  at  87  Knots  TAS) 
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Figure  38.  Response  to  Me  Mured  Airloads  (13  Modes) 
Condition  17  (Autorotation  at  83  Knots  TAS) 


Figure  k2.  Response  to  Measured  Airloads  (13  Modt:3 ) , 
Condition  2k  (Level  Flight  at  12U.5  Knots  IAS) 


Figure  M.  Response  to  Measured  Airloads  (13  Modes} 
Condition  29  (Level  Flight  at  215.5  Knots  TAS) 
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Figure  5^ .  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  1  (Hover) 


ATMAN  MB'- MHOS 

Figure  55-  Normal  Bending  Response 
+  -  Measured  Airloads  (5  lioraal 
t  Bending  Modes),  Condition  N 
(Collective  Pullup  at  0  Knot  TAS) 
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Figure  56.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  5 
(Forward  Flight  at  51  Knots  TAS) 
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Figure  57.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes ) ,  Condition  8 
(Forward  Flight  at  105  Knots  TAS) 
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Figure  5o.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  11  (Left 
Turn  at  8U  Knots  TAS) 
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Figure  59.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  12  (Right 
Turn  at  82  Knots  TAS) 
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Figure  62.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  21 
(Flare  at  60  Knots  TAS) 
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Figure  63.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  23 
(Level  Flight  at  109  Knots  TAS) 
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Figure  66.  Normal  Bendirg  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes ) ,  Condition  27 
(Level  Flight  at  227  Knots  TAS) 
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Figure  67-  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  31 
(Level  Flight  at  232  Knots  TAS) 
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Figure  68.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes ) ,  Condition  33 
(Level  Flight  at  157  Knots  TAS) 
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Figure  69-  Normal  3ending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  36 
(Puliup  at  126  Knots  TAS) 
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Figure  70.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  37 
(Pullup  at  121*  Knots  TAS’ 
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Figure  71.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes ),  Condition  39 
(Pullup  at  206  Knots  TAS) 
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Figure  72.  Hemal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  ho 
(Puilup  at  2Co  Knots  TAS) 


AZMOmMCU-MHSS 

Figure  73.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  1*3 
(Puilup  at  163  Knots  TAS) 


72 


UMffl'MSU-NMBS 

Figure  76.  Normal  Bending  Response 
to  Measured  Airloads  ( 5  Normal 
Bending  Modes),  Condition  !»9 
(Right  Turn  at  l6^  Knots  rAS) 
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Figure  77.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  5C 
(Right  Turn  at  208  Knots  TAS) 
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Figure  73.  liormal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes ) ,  Condition  2  ( Collective 
Pullups  at  0  Knot  TAS) 


Figure  79.  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes } ,  Condition  3  ( Collective 
Pullups  at  0  Knot  TAS) 
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Figure  80.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  6  (Fo.-ward 
Flight  at  59.5  Knots  TAS) 


Figure  8l.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  7  (Forward 
Flight  at  80.5  Knots  TAS) 
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Figure  82.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 


Modes ) , 

Condition  9  (Left  Turn 

at  6l  Knots  TAS) 
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AZIMUTH  AMU  -  KOTOS 

Figure  84.  Normal  3ending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  13  (Collective 
Pullup  at-  84.5  Knots  TAS) 


-—TEST 

'  ttcoay 


*  I  I  I  I 

0  T0  MB  270  MG 

AZIMUTH  AMU  •  LEGATES 


Figure  83.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  10  (Right  Turn 
at  58  Knots  TAS) 
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Figure  85.  Normal  Banding  ResTvmse  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  l4  (Collective 
Pullup  at  86  Knots  TAS) 
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Figure  86.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  15  (Collective 
Pullup  at  87  Knots  TAS) 
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Figure  88.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  18  (Transition) 
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Figure  87.  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  17  (Autorotation 
at  83  Knots  TAS) 
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Figure  89.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  20  (Transition) 
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Figure  90.  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes) ,  Condition  22  (Flare) 
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Figure  92,  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  28  (Level  Flight 
at  170  Knots  'IAS) 


Figure  91.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  2 1*  (Level  Flight 
at  128.5  Knots  TAS) 
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Figure  93.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  29  (Level  Flight 
at  215.5  Knots  TAS) 
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Figure  9h.  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  30  (Level  Flight 
at  219.5  Knots  TAS) 
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Figure  96.  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  35  (Level  Flight 
at  219  Knots  TAS) 
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Figure  95-  Normal  Bending  Response  to 
Measured  Airloads  (5  Normal  Bending 
Modes),  Condition  3N  (Level  Flight 
at  202.5  Knots  TAS) 
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Figure  97.  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  38  (Pullup 
at  160  Knots  TAS) 
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Figure  98.  Normal  Bending  Response  to  Figure  99-  Normal  Bending  Response  to 
Measured  Airloads  ( 5  Normal  Bending  Measured  Airloads  ( 5  Normal  Bending 
Modes),  Condition  Ul  (Pullup  Modes),  Condition  42  (Pullup 

at  83  Knots  TAS)  at  166  Knots  TAS) 
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Figure  100.  Normal  Bending  Response  Figuie  101.  Normal  Bending  Response 
to  Measured  Airloads  (5  Normal  to  Measured  Airloads  (5  Normal 

Bending  Modes),  Condition  U5  Bending  Modes),  Condition  Vf 

(Left  Turn  at  124  Knots  TAS)  (Left  Turn  at  207-5  Knots  TAS) 
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Figure  .1.02.  Normal  Bending  Response  to  Measured  Airloads  (5  Normal 
Bending  Modes),  Condition  43  (Right  Turn  at  124  Knots  TAS) 
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Figure  103.  Measured  and  Computed  Figure  104.  Measured  and  Computed 
Airloads,  Lockheed  Program,  Airloads,  Lockheed  Program,  Condition 

Condition  1  (Hover)  4  (Collective  Pullup  at  0  Knot  TAS) 
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Figure  111.  Measured  and  Ccr.puted 
Airloads,  Lockheed  Program, 
Condition  23  (Level  Flight 
at  109  Knots  TAS) 


Figure  112.  Measured  and  Computed 
Airloads,  Lockheed  Program, 
Condition  25  (Level  Flight 
at  163.5  Knots  TAS) 
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Figure  113.  Measured  and  Computed 
Airloads,  Lockheed  Program, 
Condition  26  (Level  Flight 
at  207  Knots  TAS) 
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Figure  11**.  Measured  ai*d  Computed 
Airloads,  Lockheed  Program, 
Condition  27  (Level  Flight 
at  227  Knots  TAS'' 


Figure  115*  Measured  and  Computed 
Airloads ,  Lockheed  Program , 
Condition  31  (Level  Flight 
at  232  Knots  TAS) 


Figure  116.  Measured  and  Computed 
A: rloads ,  Lockheed  Program, 
Condition  33  (Level  Flight 
at  157  Knots  TAS) 
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Figure  125.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  5  (Forward  Flight 
at  51  Knots  TASy 


Figure  126.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  8  (Forward  Flight 
at  105  Knots  TAS) 
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Figive  127.  Measured  and  Computed  F.gure  128.  Measured  and  Computed 
Airloads,  Cornell  Program,  Airloads,  Cornell  Program, 

Condition  11  (Left  Turn  Condition  16  (Autorotation) 

at  8L  Knots  TAS) 
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Figure  135.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  31  (Level  Flight 
at  232  Knots  TAS) 
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Figure  136.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  •  Level  Flight 
at  157  n  IAS) 


Figure  137.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  36  (Pullup 
at  126  Knots  TAS) 


Figure  133.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  37  (Pullup 
at  12h  Knots  TAS) 
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re  139-  Measured  a  :.d  Computed 
Airloads,  Cornell  Program, 
Condition  39  (Fillup 
at  206  Knots  TAS) 


Figure  140.  Measured  and  Commuted 
Airloads,  Cornell  Program, 
Condition  ^0  (Pullup 
at  206  Knots  TAS) 
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Figure  ll*5.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  6  (Forward  Flight  at 
59.5  Knots  TAS) 


Figure  ll*6 .  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  7  (Forward  Flight  at 
80.5  Knots  TAS) 
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Figure  11*7.  Measured  and  Computed 
Ai 'loads,  Cornell  Program, 
Condition  9  (Left  Turn  at  6l 
Knots  TAS) 
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s  lU8.  Measured  and  Computed 
Lrloads,  Cornell  Program, 

3  It  ion  J.0  (Right  Turn  at  58 
Knots  TAS) 
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Figure  li»9-  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  12  (Right  Turn  at  82 
Knots  TAS ) 


Figure  150.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  13  (Collective  Puliup 
at  8^.5  Knots  TAS) 
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Figure  155.  Me&surea  and  Computed  Figure  156.  Measured  and  Computed 
Airloads,  Cornell  Program,  Airloads,  Cornell  Program, 

Condition  20  (Transition)  Condition  22  (Flare) 
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Figure  157-  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  2b  (Level  Flight  at 
12b. S  Knots  TAS) 


Figure  158.  Measured  on d  Computed 
Airloads,  Cornell  Program, 
Condition  28  (Level  Flight  at 
170  Knots  TAS) 
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Figure  159.  Measured  and  Computed  Figure  l60.  Measured  and  Computed 
Airloads,  ComeU  Program,  Airloads,  Cornell  Program, 

Condition  29  (Level  Flight  at  Condition  3Q  (Level  Flight  at 

215.5  Knots  IAS)  219.5  Knots  TAS) 


Figure  l6l.  Measured  and  Computed  Figure  162.  Measured  and  Computed 
Airloads,  Cornell  Program,  Airloads,  Cornell  Program, 

Condition  31*  (Level  Flight  at  Condition  35  (Level  Flight  at 

£02.5  Knots  TAS)  219  Knots  TAS) 


Figure  163.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  38  (Pullup  at 
160  Knots  TAS) 


Figure  164.  Measured  and  Computeo 
Airloads,  Cornell  Program, 
Condition  4l  (Pullup  at 
83  Knots  TAS) 
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Figure  169.  Measured  and  Coaputed 
Airloads,  Cornell  Program, 
Condition  1*7  (Left  Turn  at 
207.5  Knots  TAS) 


Figure  170.  Measured  and  Computed 
Airloads,  Cornell  Program, 
Condition  U8  (Bight  Turn  at 
121*  Knots  TAS) 


“igure  173.  Airloads  Computed  from  Figure  17U.  Airloads  Computed  from 
Equilibrium,  Condition  8  Equilibrium,  Condition  25 

(Forward  Flight  at  105  Knots  TA3)  ( L'.-vel  Flight  at  163-5  Knots  TAS) 
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Figure  175.  Airloads  Computed  from  Figure  176.  Airloads  Computed  from 
Equilibrium,  Condition  31  Equilibrium,  Condition  37 

(Level  Flight  at  232  Knots  TAS)  (Puilup  at  12U  Knots  "AS) 
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III  oi'  a  3-volume  report 
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Fort  Eustis,  Virginia 

This  report  presents  the  results  of  a  two-phase  research  program  consisting  of 
(1)  In-flight  me  asurement  of  aerodynamic  pressures  and  structural  loads  on  a 
compound,  rigid-rotor  helicopter, and  (2)  correlation  of  these  data  with  theoretical 
results.  y 


Flight  test  data  obtained  .n  Phase  I  and  recorded  on  an  oscillograph  were  read  on  an 
oscillograph' reading  machine  and  were  processed  in  an  automatic  data  reduction 
program.  '  Th*s  d.ita  processing  consisted  of  integration  of  the  pressure  data  to 
obtain  the  distribution  of  aerodynamic  lift  and  pitching  moments  over  the  rotor 
blade,  as  functions  of  azimuth  position.  Airload  and  structural  load  data  were 
harmonically  analyzed. 

Output  of  the  data  reduction  program  was  used  in  Phase  II  as  input  to  the  correlation 
program.  The  measured  airloads  were  used  to  compute  the  theoretical  bending 
and  torsion  responses  of  the  blade.  The  measured  torsion  moments  were  used 
in  the  theoretical  prediction  of  the  airloads. )  The  results  of  the  applied  theories 
are  compared  with  the  flight  measurements.'^ 


